Preprint: ADP-97-13/T250 
riep-ph/ 97050211 



Optimization of Monte-Carlo calculations of the 
effective potential 



"|> ! A. Ardekani a ' fe * and A. G. Williams"' 6 * 

. "Department of Physics and Mathematical Physics 

b Special Research Center for the Subatomic Structure of Matter, 
c3 ■ University of Adelaide, SA 5005, Australia 

Abstract 

We study Monte Carlo calculations of the effective potential for a scalar 
q . field theory using three techniques. One of these is a new method proposed 

and tested for the first time. In each case we extract the renormalised 
quantities of the theory. The system studied in our calculations is a one 
component 4 model in two dimensions. We apply these methods to both 
the weak and strong coupling regime. In the weak coupling regime we com- 
Ph! pare our results for the renormalised quantities with those obtained from 

two-loop lattice perturbation theory. Our results are verified in the strong 
coupling regime through comparison with the strong coupling expansion. 
^ ■ We conclude that effective potential methods, when suitably chosen, can 

be accurate tools in calculations of the renormalised parameters of scalar 
field theories. 



* aardekan@physics . adelaide . edu . au 

* awilliamOphysics . adelaide . edu. au 



1 



I. INTRODUCTION 



An understanding of the underlying vacuum structure of a quantum field theory is essential 
for understanding its physical content. This analysis is conveniently carried out by calculating 
a quantity known as the effective potential denoted by U((f>) and the minimum of which 

gives information as to the nature of the lowest energy eigenstate of the theory. This makes 
U (0) very useful, particularly in studies of spontaneous symmetry breaking (SSB). The effective 
potential determines the one particle irreducible (1PI) vertices |]] at zero momenta and reflects 
any non-trivial dynamics. It is also widely used to study radiative corrections in quantum field 
theories || . Truncating the loop expansion of the effective potential often gives it a complex and 
nonconvex character, in spite of the fact that on general grounds the effective potential must be 
real and of convex character 0] . It has been pointed out that the loop expansion for the effective 
potential fails for the fields in just those regions where the classical potential is nonconvex; the 
most familiar case corresponds to a double- well potential |J. Therefore, it is important to 
carry out nonperturbative studies which can be used even where the loop expansion is not 
applicable. One convenient nonperturbative approach is to employ a discrete version of the 
theory, i.e, lattice field theory. Lattice field theories have an ultraviolet (UV) regulator (the 
lattice spacing) and and an infrared (IR) cut-off (the lattice size) and are conveniently studied 
using Monte Carlo (MC) methods. 

The model used in our study is the A0f +1 model. The Higgs mechanism is based on a 
more elaborate version of such a model and is usually discussed at the tree level. A fully 
nonperturbative treatment of the Higgs model would be of considerable interest, but is not 
discussed further here. 

In this work we investigate the effective potential for \<p\ , x using three different methods: 
The variation of the source method (VSM) || and two types of constraint effective potential, 
which we denote as CEPI and CEPII, (see refs. and @ respectively). 

We point out the advantages and disadvantages of each method and their accuracy. Some 
suggestions for improving the methods are also put to the test. In Sec. 2 we briefly summarise 
the model to be studied. In Sec. 3 we review the above methods of calculation of the effective 
potential. In Sec. 4 we perform the calculations for both the symmetric and the spontaneous 
symmetry breaking cases in the weak coupling regime and we compare our results with those 
obtained from lattice perturbation theory. We also perform calculations in the strong coupling 
regime and compare these with the strong coupling expansion. 
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II. THE A0 4 MODEL 



We start with the action of a single component A0 4 theory in d- dimensions in Euclidean 
space in the presence of a source J, (in units where H — c — 1) 



S{<P, J] — J d d x l -{d^f + im 2 2 + ±<f _ J(j) 



A discrete lattice version of the action can be written as 



S[$, J] 



\ E(^,M " + \ E + E T^n " E ^ 



(1) 



where we have defined the dimensionless quantities = a^/ 2 ^ _1 0, m = ma and A = \a 4 ~ d 
and J = a( d / 2 ) +1 J. In addition n = (n 1: . . . ,n d ) is a d-dimensional vector labeling the lattice 
sites and ji is a unit vector in the temporal or spatial direction. The sum over jj, is over the 
d Euclidean directions. We also have denoted the field on the neighboring site of n in the 
direction of fi by <fi n ,n- Henceforth we drop the hat from the dimensionless field variables and 
sources for brevity unless it is necessary to avoid confusion. We also impose the appropriate 
periodic boundary condition on fields: 



for all /I , (2) 



where JV M = (0, . . . , . . . , 0), is a d-dimensional vector with being the number of lattice 
sites in the direction \i. 

The 4 theory is known to exist in two phases, one where the reflection symmetry <p — > — <fi 
is spontaneously broken and the other where it is not. The symmetric phase with (0) = 0, 
is separated from the broken symmetry phase with ((f)) ^ by a line of second order phase 
transitions where m and A assume the critical values m c and A c . 

For the action S[<f>] on the lattice the generating functional for the correlation functions is 
defined as: 

J[d0]e- 5[< ^ J] 

Z[J] = JM^ ] ' (3) 

such that Z[0] = 1. /,From Z[J] one can define the connected Green's functions as: 

G(m, . . . , nj ) c = _..._ w[J}\ J=0 , (4) 

""Til n j 

where 

W[J]=]nZ[J]. (5) 
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III. THE LATTICE EFFECTIVE POTENTIAL 



Consider a lattice lagrangian density on a <i-dimensional cubic lattice with the total number 
of lattice sites N d , 

£n = Y.\^-<Pnf + V(<P n ) . (6) 

The classical vacuum ( ground state) is at the minimum of V{4>). The vacuum expectation 
value {(f)) of the quantum field is not necessarily identical to the classical vacuum. The vacuum 
expectation value of the field in the presence of an external source, J(x) is given by 

The vacuum expectation value (0) is the limit of C „ as J — > 0. Hence we can ask for what 
value of J can one obtain a given C . One can choose to treat C as the independent variable 
instead of J and define the "effective action" T[0 C ] by a Legendre transformation: 

T[(j )c }=J2^nJn-W[J}, (8) 

n 

where C is defined in Eq. (|7|). It is easy to verify that 

In the case J = 0, by translational invariance it follows that C must become constant (i.e, 
independent of the label n). Hence the vacuum expectation value is given by (0) and satisfies 

dT\<p c ] 



= . (10) 

--(4>) 

Similarly for any constant J we must have C = <fi also constant. Define the effective potential, 
U@), by 

r$ = N d ug>) . (ii) 



The Fourier trasform on a finite, discrete lattice is defined by 

fe ^^e 2 ™^0 n , (12) 

n 

where k = hi, . . . , kd is a d dimensional vector with {—N/2) < k n < N/2, (we assume N is even 
from this point) and where n.k = ri\ki + . . . + ndkd. The coordinate-space and momentum-space 
5-functions are 
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X — 1 -2m(n-m).k/N r J_ -2m(k-q).n/N Qo\ 



fc 

respectively. The inverse Fourier transform is 

Note that we have used the asymmetric normalisation of the Fourier transform and its inverse 
as is usual in the field theory in the continuum. The effective action is the generator of proper 
(i.e, one-particle irreducible) Green's functions and in particular we can Taylor expand the 
effective action to give 

OO 1 

n^HETn E rW(n 1; ...,n M )0 cni ...0 cnM . (14) 

m=o M ■ m,...,n M 

Here T^ M \ni, . . . ,Um) is the proper M-point Green's functions in presence of the source J n 

d»m =r( Ao (ni> ... )njtf) , (15) 



'~"rcn\ ■ ■ ■ '-'Yan.M 

In terms of its Fourier transforms we have 

00 1 1 ~ ^ ^ ~ ~ 

m)=E w ^ d E fw (fcl) ...M^-^ M » as) 

M=0 ' k u ...,k M 

where here is the Fourier transform of C 

The vacuum proper Green's functions are obtained by setting J = 0. If the source is a 
constant (J n = J for all n), then the translational invariance is restored and we can factor out 
an overall normalisation and a ^-function to define 

f ...,k u ) = N\^ i+ _^Ji M Hh, k M ) (17) 

Where in the limit J — ► we recognise that f( M ) (hi,..., k^) is the dimensionless, lattice 
equivalent of the proper M-point Green's function in momentum space. For constant J we also 
have 4> cn — > <fi and hence = N d 5^ Q and Eq. (0) gives 

m = N d Y: ^fi M) (o)f, (is) 

where 

fW(0)=fW(0,...,0). (19) 
Comparing Eq. (|T~TJ) with Eq. ([16]) gives 
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It immediately follows that 



W) = E t^ m) W 



d M U{(t)) 



—rM 



(20) 



T M 



= ff)(0), (21) 

<M$=0 

where here it is understood that we are working in the unbroken symmetry sector, (0) = 0. In 
the unbroken sector we see from Eq. (|20"D that the dimensionless, proper Greens functions with 
vanishing momenta can be easily obtained from the effective potential, U(4>), by differentiation. 
We see that minimises U(<p) and in the limit J — > the minimum — > (0). Also note that 
Eq. (|[) gives an expansion of J in terms of the 0's and r(0)'s 

1 



-fW(0) 



22) 



m=AM-1)\ 

In the broken symmetry sector, (0) ^ 0, it is more appropriate to use the shifted field 

X (x) = <t>{x) - (0) . (23) 

The one-particle irreducible (IPI) vertex functions are linear combination of the r^ M ^'s, 
and can be obtained from the shifted version of Eq. ( |2T)| ) 

x ^-(0}] M 



f/(0) = f/ (s) (x) = f/ (S )(0-(0))= E ■ 

M=0 



M! 



(0). 



(24) 



As is usually done in lattice field theory studies we renormalise at the renormalisation point 
where all external momenta of the Greens functions vanish. The renormalised quantities can 
be obtained directly from the effective potential. For example in the A0 4 theory we have: 

dU($) 







d 2 (f) 



72 d 4 U(<p) 



zrf (o) = rf 



<4>) 



z 2 r^(o) = r^ = A r 



(25) 



(26) 



(27) 



where y Z is the field wavefunction renormalisation constant (0 r = y Z(p). ^From the first two 
conditions above and requiring m 2 . > it follows that (0) is at the minimum of U(<fr). Also note 
that m r and A r defined as above are not the physical mass and coupling, which are defined in 
the pole of the propagator in the complex energy plane and the on shell four-point function, 
respectively. However in the scaling region (close to the critical line) these values are a good 
approximation to the physical mass and coupling M. 
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IV. THE MC EFFECTIVE POTENTIAL. 



In this section we will examine three MC methods for calculation of the lattice effective 
potential. The renomalised coupling constants obtained by these methods are compared with 
analytical results. From this point on we work exclusively in two dimensions (d=2). 

A. The variation of source method (VSM) 

Eq. (^2|) suggests that in the Monte Carlo calculation one can calculate the mean value of 
the fields, 0, for different values of the source and as a result one obtains as a function of J. 
This function can then be inverted to obtain J as a function of 0, i.e., J(0). Then using Eq. ( |22|) 
we see that the derivatives of J with respect to would give the proper Green's functions at 
zero momentum. ^From Eq. ([/]) one also concludes that 0j is antisymmetric in J. That is 

4>j = • (28) 

Fig. 1 shows J(4>) as a function of for the symmetric case (Fig. la) and the broken symmetry 
case (Fig. lb). Note that for the broken symmetry case, 0j as a function of J is discontinuous 
and so the relation in Eq. (p2|) can not be inverted for all 0(J). Whenever it is possible Eq. (|22|) 
has to be inverted to obtain the source J as a function of 0. Then the derivatives of J with 
respect to would give the vertex functions at zero momenta and consequently the renormalised 
masses and couplings can be calculated. 

The mean value of the field in the presence of a source has a small statistical error. This is 
expected since it is an analog to the reduction of fluctuations of a spin system in the presence 
of an external magnetic field. As the source becomes smaller the fluctuations become larger. 
Thus one needs to perform the calculations for large enough sources that the error is small and 
then extrapolate the results to J = 0. 

This method will be referred to as the variation of source method (VSM) and has a number 
of advantages. The vacuum expectation values of the field 0(J) are the simplest quantities 
to compute on the lattice and their J-dependence can be exploited to get the first derivative 
of the effective potential. Since the source effectively causes the boson field to become more 
massive, the finite size effects generated by the lattice become exponentially small provided 
that the lattice is large enough. Since the data become noisy for small values of J we need 
to restrict the analysis to a safe region of J, which can introduce some errors in the results 
through uncertainties in the extrapolation. 
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B. The Constraint Effective Potential (CEP I) 



In the last section the effective action and the effective potential U (0) were defined through 
introduction of a source J. There is a different method which does not require such a dy- 
namical symmetry breaking source. The constraint effective potential was first introduced by 
Fukuda and Kyriakoloulos J7j as an alternative way of obtaining the explicit expression for 
the effective potential. It was further analyzed by O'Raifertaigh, Wipf and Yoneyama 0. In 
this approach one obtains an explicit expression for the effective potential, without introducing 
external sources, but instead through the introduction of a ^-function in the functional integral. 
In the constraint effective potential approach one first defines U(<p) = U(N 2 ,(f)) as 

e~ N2m = J [d<P}6 (± £ 0n - ^ e-W (29) 

and then uses the fact that as iV 2 — > oo we have U(<f>) — > U{<f>) and the effective potential is 
recovered. W 

It is easiest to demonstrate this result in Minkowski space, where Eq. ( |29|) becomes 

e 

We can replace the 5-function in Eq. (|3~0"D by its integral representation to obtain (up to an 
irrelevant constant) 

e -i*m = c J dj J [ d(f) ] e ii^+J4-]-iN^j 

= C fdJ e^- N2j ^ (31) 



" X ' ll[ ~ ] = '[W^E^-^M. (30) 



Note that in the integrand of Eq. (plf) we have (j) fixed and J arbitrary. In the limit, N — > oo 
the dominate contribution to the integral comes from the stationary point of the integral which 
is the value of J at which dW[J]/dJ = 0. Recall that T(0) = (J0- W[J}) 
which we see that up to an irrelevant overall constant 



, from 

(j,=dW[J]/dJ : 



e -iN>U(<j>) _ e -iT W = e -iN*U W ag ^2 ^ QO) (32) 

as claimed. 

We can also arrive at this result directly in Euclidean space by multiplying both sides of 
Eq. ( p9|) by e N2j< ^ with J arbitrary and then integrating over <p to obtain 

fy e -N*m-Jft = f[d<p} e - s ^ +J ^^. (33) 



As iV 2 — > oo the left hand side of Eq. ( p3|) becomes entirely dominated by the stationary point 
of the one- dimensional integration given by dU(<f))/d<f) = J, while the right hand side is 
recognised as e w ^ for a constant source, J. Hence up to an irrelevant overall constant we find 



8 



e -N*M4>)-J4i _ e W[J] as ^2 ^ QO) (34) 
and so find that (up to a constant) 

e -**U® ^ e W[J\-N*J$ = e -N*U® as ^2 ^ QO) (35) 

as required. 

It is important to note that the e~ N2u<y( ^ relates to similar definitions in statistical mechanics 
and spin systems [10] and that 

e -N 2 U($) 

P ^ ~ Jd~4>e- N2 uil>) ' (36) 

can be interpreted as the probability density for the system to be in a state of "magnetization" , 
<p. Then it can be seen that the probability for a state whose average field is not a minimum of 
U(<f) then decreases as N 2 — ► oo. 

This suggests that one needs to study the probability distribution of the order parameter 
(f). Using a Monte Carlo algorithm one generates a Boltzman ensemble of configurations, {</>}, 
weighted by e~ s ™. Let dj\f be the number of configurations with average field values in an 
interval dip about 0. Then 

dM{4>) = Ce- N ^d4> , (37) 
with C some constant. Then one can write 

up to an irrelevant additive constant. Eq. (|36|) suggests that one can generate a large number of 
configurations weighted by e~ 5 ^ , calculate <fi for each configuration and construct a normalised 
histogram. The histogram can be fitted to Eq. (|38D . The most probable average field values 
are near the minimum of the effective potential. In order to determine U(<p) away from its 
minimum, i.e., to sample a range of relatively improbable values of 0, one can introduce a small 



source. Then a simple generalisation of Eq. (|36|) allows a nonzero external source |TT 



Thus one can check whether such an ansatz gives a good approximation for the effective po- 
tential, and so construct the effective potential by performing a simultaneous fit of several 
histograms corresponding to different values of J. By the expression "simultaneous fit", we 
mean that the chi-squared values corresponding to each J are summed and this sum is then 
minimised. This method can be applied easily on the lattice. Note that in Eq. (^) we have as- 
sumed that for sufficiently large iV 2 the finite volume effects on U{4>) can be neglected, i.e., that 
the lattice volume is sufficiently large. The constraint effective potential method summarised 
in Eq. (|39|) will be refered to as CEPI. 
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C. The Constraint Effective Potential (CEP II) 



Now return to Eq. ( p8|) and perform a shift of field, 4>(x) — > 0(x) + 0. Since the measure is 
translationally invariant we obtain 



e 



- ^M(^?£0») e " 5[ ^ (40) 



Taking the derivative with respect to <fi we get 



N>m = 1 / [#] ( y ( % ) ^ + ^) c -*[*3 (41) 



# 6 V 2 7 L ^ J " v A^ 2 

Only the potential part of the action is affected by the shift of field since is constant and so 
dS/dcf) = N 2 dV/d(p. Using this fact and shifting the field back to its original form then gives 



dU(J>) _l dV{<f>) \ 
dcj) \ d<p /- 



(42) 



where we have introduced the shorthand notation 



{Om-$=i^ m ) J Wi 5 (^£^ -<« 0(<f>)e- s M. (43) 
In the A0 4 theory being considered here we find 

^M=m 2 + Vh- (44) 

Expressions for some of the higher derivatives of U(4>) are given in Appendix A. These equations 
are very useful in the Monte Carlo (MC) calculations since they relate the derivatives of the 
effective potential (and consequently the zero momentum vertex function) to the averages of 
quantities that can be calculated directly from the lattice. This method will be refered to as 
CEPII. 

There are two ways of calculating the renormalised quantities using CEP II. The first one 
applies the constraint on the lattice, fixing 0, then calculates (0 3 ), and finally uses Eq. (|44"D 
to obtain the first derivative of the effective potential. Higher derivative are evaluated from 
fitting a curve to the dU/d<j) versus results. This has some similarities with the variation of 
source method, however there is a difference between these two methods. In VSM one sets the 
source J to constant and (4>) = up to fluctuations due to finite A^ 2 , whereas in CEPII we 
have ((p) = (J) exactly by construction. 

In the broken symmetry sector there is another difference between this method and VSM 
in the broken sector. When using VSM we are not able to obtain any value of (0) in the region 
between the two minima, whereas both CEP methods are suitable for probing this region. One 
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can always fix <\> to any value including the values between the two minima to get the full shape 
of J{4>) (see Fig. 2). However as far as the practical calculation of renormalised quantities is 
concerned, this method is almost equivalent to VSM and so from here on we disregard this 
approach. 

The second approach to CEPII is through the equations shown in Appendix A and is more 
direct. These equations relate the derivative of the effective potential to the averages of some 
derivatives of the classical potential. All these averages should be taken in the presence of the 
constraint which fixes {(f)) = (p. 

Imposing a constraint on a Monte Carlo algorithm is relatively simple. One starts with a 
configuration with the field average being <fi [0. Each time a site is updated by a value 5 (say) 
such that 

<f>i = & + 5 , (45) 

then some arbitrary site k must be updated simultaneously such that <fi' k = tpf. — 5. This 
procedure is carried out for all the sites, which completes a sweep. For a large enough lattice 
imposing such a constraint is not expected to violate the ergodicity of the MC algorithm. 

The advantage of this method over the VSM is that one does not need to run a Monte Carlo 
routine several times with different sources, and no curve fitting is required. One disadvantage 
of this method is that for calculation of the renormalised coupling one needs to add and subtract 
many average terms as has been shown in Appendix A. Although the statistical errors might be 
small for each term, the overall errors contributing to the renormalised coupling can be large. 
However, the renormalised mass in the symmetric phase of the A0 4 theory obtained using this 
method is very accurate. 

We also would like to comment on Fig. 2. It has been shown by a very general argument 
that U"(4>) > for all <fi [|J, (primes denote differentiation with respect to 0). This general 
property is known as the " convexity" of the effective potential. Looking at Fig. 2 it is clear 
that this condition is violated for U(<f)). This can be understood by noting that convexity holds 
only in the thermodynamical limit, i.e., N 2 — > oo. 

To conclude this section it should also be mentioned that the proper vertex functions can 
be obtained directly using the standard Monte Carlo (MC) method. For example for the A0 4 
four point vertex function one obtains 




Here, for example, (0 4 ) c is the connected part of vacuum expectation value of fourth power of 
the Fourier transform of the field at zero momentum. As we will show, in the weak coupling 
regime this method suffers from very noisy signals giving rise to large statistical errors. The 
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errors are due to the large fluctuations of correlation functions in this regime as well as the 
subtraction of the disconnected pieces. However in the strong coupling regime this method gives 
a relatively good approximation for Ty)(0) and the statistical errors are reasonably small [12 



However the higher order vertex functions calculated with this approach can be very noisy even 
in the strong coupling regime, primarily due to subtractions of noisy disconnected pieces. 



V. THE NUMERICAL RESULTS 

In this section we present our results for the calculation of renormalised coupling, X r , in two 
dimensions. It includes the symmetric and broken symmetry sector in the weak coupling regime 
as well as the strong coupling regime. In the case of the weak coupling regime the results are 
compared with 2-loop results and the direct calculation of A r using the standard MC method 
in Eq. (f46|). In the strong coupling regime we also compared the results of each method with 
the strong coupling expansion results. The details of the numerical simulation are included at 
the end of this article. 



Case 1: The symmetric sector in the weak coupling regime (WCR) 

The variation of the source method (VSM): 
Here we study the model in the symmetric sector where ((f)) = 0. As we will see, all methods 
presented in this paper require the calculation of renormalised mass m r , and the wavefunction 
renormalisation constant, Z. In general, the boson propagator extracted from the lattice has 
the form 

G(k) = ^ |— ^- , (47) 

k 2 + m 2 {k 2 ) 

where m r = m r (rn^) is the mass-pole of the scalar particle, i.e., the renomalised mass. In 
particular at zero momentum 

G(k = 0) = ^ , (48) 

where we make the standard approximation that m 2 « m 2 (0). The renormalised mass, m r , is 
then given as the reciprocal of (, the lattice correlation length 



^2 _ 1 1 dG{k) 

G(k) dk 2 



(49) 

fc=0 
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Taking into account the translational invariance of the correlation functions, one can choose to 
approximate the momentum derivative in Eq. fl49|) by the variation of G(k) across one lattice 
spacing and in one direction to obtain 



N 
2^ 



( 7 2 ) c - (a 2 ) c - ((3 2 ) c 
(« 2 )c + (/9 2 ) c 



(50) 



with 



N N 



a 



EE 

n=l m=l 
N N 

EE 

n=l m=l 



J n,m 



COS 



sin 



27T 

— (n 

N 2' 



,2tt N 
— [n 



7 



N N 

EE 

71=1 771=1 



J n,mi 



(51) 



where here n, m label the temporal and spatial coordinates for the field <fi respectively. 

There are two different ways of calculating Z. One is to use Eq. and the fact that 
G(0) = N 2 (<ft 2 ) to calculate Z. The second way of calculating Z comes from combining Eqs. 
and (]27|) which gives 



d 2 U{(t)) 



G(0) 



(52) 



Thus G(0) can be directly calculated from the fit and the calculation of Z follows as before. 
An accurate calculation of m r is crucial for both methods. We found that in the weak coupling 
regime, the second method was more precise. We compared our results with the 2-loop lattice 
perturbation theory calculations (LPT) of the renormalised parameters. This means that finite 
size effects may be present in our comparisons at some level. The comparison is shown in Fig. 
3. The values for A r seem to be accurate even in the very weak coupling regime. In this regime 
the effective potential results are in good agreement with the lattice perturbation calculations. 
The MC results begin to deviate from the perturbative calculations as A r increases. This is 
expected since a loop expansion in A0 4 theory is an expansion in A r and as this is increased the 
contribution from higher loops becomes more significant. 

The VMS can be expensive in CPU time but the cost can be reduced to some extent. For 
a value of J it is possible to calculate D n (<fr) = d(j)j/d n J = N 2 ((j) n ) c during the calculation of 
(fij, for each value of J. From these derivatives one can expand <pj around J and then use a 
curve fitting routine to calculate A r , as we did before. The statistical errors can become larger 
for the higher derivatives because of the subtraction of the disconnected pieces of D n (<ft). In 
Table 1 we have shown a comparison of our previous results for 0j and results obtained by 
expansion around three J values, namely, J = 0.075,0.25,0.4 for m 2 = 0.1, A = 0.1. We see 
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that the calculated values of <p are reasonably close to the previous results. However, the price 
for reduced computational time is a slight increase in uncertainties. 

We have also calculated A r for A = 0.055, m 2 = 0.1 using Eq. (fHf ) and the result is included 
in Fig. 3. The statistical errors are extremely large and it suggests that the calculation of the 
4-point vertex function in this region is impractical with this method. 

The constraint effective potential method I (CEPI): 

This method is the easiest to implement. We generated the Boltzman ensemble of indepen- 
dent configurations. For every configuration we measured <ft — (l/-^ 2 ) an d computed the 
histograms for the probability density P(<p) for several values of J. We also noticed that the 
anzatz of Eq. (|39|) only worked well for very small J in this region. We did a simultaneous fit 
to Eq. fl3"6|) of a few histograms corresponding to J = and small J's using a three-parameter 
anzatz for U((f>) of the form 

U (4>) = ai0 2 + a 2 4 + a 3 6 . (53) 

Although there was no systematic discrepancy between the data and the fit, the statistical 
errors were very large. We unsuccessfully tried more histograms and higher powers of <fr i n the 
fit. The statistical errors remained large and we concluded that even a reasonable estimate of 
renormalised parameters in this region was not feasible with this method. 
The constraint effective potential method II (CEPII): 

In this method Eqs. ( |X4]) and (|A8|) can be used for calculations of m r and A r respectively. 
All averages shown in these equations are to be taken with the constraint of = 0. Although 
the statistical errors for each term are small, the overall error can be large. However in the 
symmetric case in in the weak coupling regime most of the terms either vanish at <fi = or are 
small enough to be neglected. For example for m r only three terms need to be considered. But 
the computation of A r suffers from larger cumulative errors. 

The results are compared with the VSM results and are shown in Table 2. We also compared 
the calculation of renormalised mass using Eq. Q5U| ) with the CEPII calculations in Table 3. The 
comparison indicates that in this sector the CEPII method can provide an accurate calculation 
of the renormalised vertex functions. 

Case 2: The broken symmetry sector 

In this section we consider the calculation of the renormalised mass and renormalised cou- 
pling in the broken sector, (0) ^ 0, in the weak coupling regime. The VSM procedure is exactly 
the same as for the symmetric sector. For fixed m 2 = —0.1 and < A < 0.17 we calculated 
m r and A r for different values of A. The error on G(0) is larger than the symmetric case due 
to the subtraction of the disconnected pieces. Thus we used Eq. ( |52] ) to calculate G(0) and 
subsequently extracted Z as previously discussed. 
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In order to calculate the renormalised quantities using lattice perturbation theory we fol- 
lowed the standard approach to treating the broken sector. That is, in the bare Lagrangian we 
shifted the field by the classical value of the vacuum, v = \J—6m 2 /\, such that 

X (x) = <f>{x) + v . (54) 

After this translation the mean value of the field, (x), vanishes and the perturbative calculation 
proceeds in the standard manner, keeping in mind that a non-symmetric x 3 interaction has been 
generated. In lattice perturbation theory one then needs to also consider vertex functions with 
a three point interaction. 

The comparison between the two-loop results and the results from the VSM method are 
shown in Fig. 4. In applying the CEPII method to the broken symmetry sector, evaluation of 
all the terms in Eq. (app8) is necessary. This renders this method impractical. As one might 
expect from the symmetric sector results the calculation of the renormalised parameters using 
CEPI also suffers from large noise difficulties and the signal could not be recovered. 



Case 3 '.Strong coupling regime 



In a weak coupling expansion the interactive term is pulled out of the path integral repre- 
sentation of the partition function as a functional operator. That is 



Z[J] = exp 



A v — 



[d<f>] exp 



E 

n,fx 



1 



V/J 2 + 7^ 2 4>n + J n(, 



(55) 



The remaining functional is Gaussian and can be done exactly. The partition function can then 
be written in terms of a power series of A and the standard perturbation theory follows. 

The strong coupling expansion was first proposed by the authors of Ref. ||13| . For this 
expansion, unlike the weak coupling expansion, the kinetic and the mass terms are pulled out 
of the path integral as a functional operator. That is 



Z[J] = exp 



m,n 5J n fiJm 



Z [J] 



(56) 



where 



Z [J] 



exp 



jfo + Jnl 



(57) 



The remaining functional integral is not Gaussian but can be evaluated as a product of ordinary 
functions on the lattice, 



F(x) 



n F(0) 



(5* 
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where 



Fix) = / dze 



J,Z 4 + XZ 



(59) 



and J\f is a constant. The function is a transcendental function and can be expanded as 
a power series in x 



j . 1 ~ 2 n a; 2n „,w 1, 

F(l) = 7cj£(*or r( 2 + 4>- 



(60) 



Using this series expansion one can easily expand the both terms in the r.h.s. of Eq. (|56|) to 
obtain a power series expansion for Z[J] which assumes the general form 



Z[J] = M' 



k=l 



(61) 



where Af-[J] are integrals over the source function J. Thus the strong coupling expansion is an 
expansion in powers of \~ k / 2 . Bender et al. fT3|J14|| obtained a series expansion for G = X r / 
of the form 



4-d 



where 



For fixed x one has 



L N 



n. .1 



G = y~ d/2 E E a m *-V 

1=0 n=0 



(62) 



x = y 



and y = ? = ± r . 



m r 



(63) 



where 



G = y-^Y, 4 

1=0 



{N) <x) y l 



(64) 



a\ N) (x) = E a inX n 



N 



(65) 



n=0 



This series does not converge for large correlation lengths. Thus the authors of Ref. 
proposed a scheme to extrapolate the expression for A r to large y assuming that A r remains 
finite in the limit y — ► 0. 

Raising Eq. fl62|) to the power of 2L/d and expanding to order L we find 



2L 
d 



y- L Y,i N \*)y l 



(66) 



\l=0 
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We then find 

/ L x (d/2L) 

G = y- d/2 (J2bi(x)y l ) , (67) 
\i=o ) 

which is equivalent to Eq. (|62|) for small y and approaches [b^ (x)] d ^ 21 in the limit y — > oo. In this 
manner the authors of Ref. |TJ| obtained an analytical series for Eq. (|67]). Since the interesting 
physics lies in a regime where the correlation length is large, we performed our calculation in 
this regime. Thus the above extrapolation scheme was necessary. 

We chose a moderate correlation length ( = 3.6 by an appropriate tuning of the bare 
parameters. This can be done by fixing A and choosing m to be in the symmetric region. As 
one decreases m, one gets closer to the critical line and the correlation length increases. Using 
this, one can reach the required correlation length. 

To apply VSM we followed the same procedure as before. For six different A's and fixed 
correlation length ( = 3.6 =F 4%, we calculated the values of <f)j for different values of J. The 
curve fitting procedure was carried out in the same way as for the previous cases. We noticed 
that in this regime the inclusion of larger 0j's can change the behavior of the fit at small <pj, the 
region which is of most interest to us. The problem arises due to the curve fitting procedure. In 
the WCR the data points close to = have much larger weighting that the one far away from 
this point. Thus, calculating the derivatives of U(<j)) at = seem to be reliable. However, 
in the strong coupling regime, the data points that are far away from = have much higher 
weighting and even a small fluctuation might affect the calculated J(0) considerably. 

We improved the results by imposing the condition in Eq. fl52|) , that is to fixing the coefficient 
ai = 1/G(0) where ai is defined in Eq. ( |53"1) and G(0) = iV 2 (0 2 ). This improved the results 
and the inclusion of larger 0/s did not affect the results significantly (up to 3%). 

Next we calculated the renormalised parameters using CEPI. Unlike the previous cases the 
errors in the results were reasonable. For the extraction of renormalised parameters we only 
used two histograms corresponding to J = and J = 0.005. In the weak coupling regime where 
the mass term is dominant, one needs to sample the higher values of i n order to improve the 
calculation of A r . Thus in the strong coupling regime there might not be a need for additional 
histograms. From the VSM results one might expect that sampling very high might have a 
similar problem. This was confirmed from our data for this particular case. 

We also calculated the renormalised coupling using Eq. (flop . Unlike the weak coupling 
regime, uncertainties in the results in this region were reasonable. All the results in the strong 
coupling regime are shown in Fig. 5 and 6. They also are compared with the strong coupling 
expansion results. They all seem to be in agreement with each other within errors. This 
indicates that as the coupling increases the MC results approach the strong coupling expansion 
results. In order to apply the CEPII method, there are numerous terms in (A. 8) which have 
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to be evaluated and consequently the accumulated errors can be very large. However we found 
that as before the renormalised mass can be calculated accurately. 

Details of the simulations 

In our MC calculation we chose the Hybrid MC algorithm. In a run we have taken a 
number of decorrelation MC iterations between two measurements. All the calculations were 
done on a 20 2 lattice and the rate of acceptance was kept between 40% and 60%. In all 
cases (except the broken sector) the calculations of renormalised mass and G(0) (where it 
was needed) and the direct calculation of T^\0) were done using 6,800 uncorrelated samples 
with 50,000 thermalisation configurations. In the broken sector we used 11,000 uncorrelated 
samples with the same thermalisation configurations. The reason for the increase was to obtain 
better statistics, since the measured quantities have larger errors due to the non-vanishing 
disconnected pieces. In applying the VSM to the symmetric case (and in the weak coupling 
regime), we calculated <fi with 0.025 < J < 0.425. We noticed that the necessary number of 
decorrelation iterations in the presence of nonzero J was smaller than for the J = case. The 
calculations were carried out using 2, 500 decorrelated configurations. We took the number of 
thermalisation configurations to be 10, 000. In the broken symmetry sector we increased the 
number of uncorrelated configurations to 3, 200. In the strong coupling regime only the range 
of the values for J was different (as mentioned in the previous section). For CEPII we used 
5, 000 uncorrelated configurations with 50, 000 thermalisation iterations. In construction of the 
probability distribution histograms, we used 750, 000 configurations. The curve fits were done 
using a standard x 2 fitting algorithm where the uncertainties on the parameters were obtained 
from the diagonal of the covariant matrix. For the strong coupling results we also estimated 
the systematic error due to the fact that ( was fixed to be approximately 4% by varying the 
fixed value within reasonable limits. 

VI. CONCLUSIONS 

We have studied the calculation of the effective potential for X(f)f +1 theory using three 
different methods: the variation of source method (VSM) and two constraint effective potential 
methods (CEP I) and (CEP II). Using our new method, referred to as CEPII, we showed how to 
calculate the vertex function using the correlation functions in the presence of a constraint. We 
calculated the effective potential in the symmetric and the broken sector in the weak coupling 
regime as well as in the symmetric sector in the strong coupling regime. The renormalised 
quantities, A r and fh r , were then obtained from the effective potential for each case. In the 
weak coupling regime we compared our results with lattice perturbation theory. We found that 
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in the symmetric case both VSM and CEP II can give accurate results, whereas the CEPI 
method and the direct Monte Carlo calculation of the (2 and 4 point) vertex functions failed 
to do so. We also found that in the broken symmetry sector VSM is the most practical and 
accurate of these methods. We also studied the model in the strong coupling regime and the 
results were compared with the strong coupling expansion results. In this regime we found that 
CEPI, VSM, and the results from the direct Monte Carlo calculation of the vertex functions 
were consistent with each other and with the strong coupling expansion results. In summary 
then, we have shown that Monte Carlo effective potential methods can be accurate and reliable 
tools for calculating physical quantities for scalar field theories, but that one should use the 
method of evaluating the effective potential and its derivatives which is best suited to the regime 
of interest. 



The differential equations relating the constraint effective potential and the classical poten- 
tial are : 

For the first derivative we have 



APPENDIX A: 




(Al) 



which for A0 4 theory becomes 



dU(<j>) 
d4> 



™ 2 + -(0 3 >. 



(A2) 



The second derivative is given by 



u(<j>) _ / d 2 v(4>i) \ 




dVjfa) dV{<j>i) \ 
#1 d<pi I 





2 



dcf \ dtf 1^ 



+ N d 



dV(fa) 
d<P 



(A3) 



which for A0 4 theory becomes 



d 2 U((j)) 

d 2 4> 




<# 3 > 



(A4) 



The third derivative gives 



19 



d"U{<t>) = Nd d 2 U_(<P) dU^) md ( dU(J>) \ 2 _ N2d f dU(j>) \ 3 / d*V X0Q 



which for A0 4 theory becomes 



dm 



(9(A 2 ) terms + 0(m 2 ) terms + 
Finally for the fourth-derivative we obtain 

?; * \m - v<i i , 

# 2 / 



terms that vanish at 6 = 



+ • 



^(0) = d (dHmV _ iN2d d 2 U(<f>) (dU(J>)\ 2 I tPVfa 
' r ' : 1 \ do ) \ do. 



d(f) 

which for A0 4 theory gives 

d 4 u(~4>) 



+ 



= A - Am 2 - iV rf Am 2 (0 2 > - N d \m 2 ^ 



+ 



0(\ 2 ) terms + 0(m 4 ) terms + terms that vanish at = 
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TABLES 



TABLE I. The comparison of the calculations of (pj for different values of J's and perturbative 
calculations of ~4>j around J = 0.1, J = 0.225 and J = 0.425 with in 2 = 0.1, A = 0.055 and N = 20. 



J 


<Pj 


error 


U 1 


error 


u.uou 


U.4UOZ 


U.UUoZ 


n a no 
u.4uy 


U.UUD 1 


U.U ( 


U.0O41 


U.UUol 


U.OooU 


U.UUOZ 


n 1 nn 
U.lUU 


n 7RA Q 
U. f 04<3 


U.UUol 


U. < t)4<5 


U.UUol 


U.lZO 


n QO^Q 

u.yzoo 


n nno7 
U.UUZ i 


U.yozU 


n nn/io 
U.UU4Z 


n i 


i ns/in 

1.U04U 


n nno7 

U.UUZ ( 


i.uyo 


n nn/i a 

U.UU48 


n 1 7^ 

U.l / D 


1 01 1 
1.Z11Z 


n nno7 
U.UUZ i 


i.zio 


n nn/i 7 
U.UU4 / 


n onn 
U.zUU 


1 QQOO 

i.ooyy 


n nno7 
U.UUZ i 


i a ^i n 
1.4D1U 


U.UUoy 


n oo^ 


i /i ^n^ 

1.4DUO 


U.UUZO 


1 /i ^n^ 

1.4DUO 


n nnofi 

U.UUZO 


0.250 


1.5590 


0.0026 


1.5646 


0.0038 


0.275 


1.6659 


0.0026 


1.6680 


0.0043 


0.300 


1.7620 


0.0026 


1.7720 


0.0049 


0.325 


1.8564 


0.0026 


1.8706 


0.0048 


0.350 


1.9411 


0.0026 


1.945 


0.0038 


0.375 


2.0166 


0.0026 


2.0209 


0.0033 


0.400 


2.0963 


0.0026 


2.0963 


0.0026 


0.425 


2.1704 


0.0026 


2.1714 


0.0029 


0.450 


2.2430 


0.0026 


2.2470 


0.0029 
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TABLE II. The comparison of the calculations of A r using the VSM and the CEPII in the sym- 
metric sector and weak coupling regime. X^ ons denotes the renormalised coupling calculated by CEPII. 
Here m 2 = 0.1 and JV = 20. 



A 


A* 


error 


\ cons 
A r 


error 


0.02 


0.0191 


0.0003 


0.018 


0.0007 




u.uoou 


u.uuuo 


\J .\JO\jd 


u.uuuo 


0.055 


0.0518 


0.0008 


0.0510 


0.0017 


0.07 


0.0670 


0.0009 


0.0657 


0.0019 


0.1 


0.0891 


0.0008 


0.092 


0.0016 


0.13 


0.112 


0.0013 


0.121 


0.004 


0.19 


0.165 


0.002 


0.175 


0.0061 


0.24 


0.216 


0.0023 


0.22 


0.007 


0.35 


0.313 


0.0035 


0.321 


0.018 


TABLE III. The 


comparison 


of the calculations of m 2 


using the VSM and the 


CEPII in the 


symmetric sector and weak coupling regime. X^° ns denotes the renormalised mass calculated by CEPII. 


Here m 2 = 0.1 and N 


= 20. 








A 


m r 


error 


m c r ons 


error 


0.02 


0.324 


0.001 


0.323 


0.007 


0.04 


0.334 


0.002 


0.330 


0.007 


0.055 


0.340 


0.0008 


0.332 


0.008 


0.07 


0.343 


0.0014 


0.339 


0.006 


0.1 


0.345 


0.0023 


0.347 


0.008 


0.13 


0.350 


0.0023 


0.357 


0.009 


0.19 


0.375 


0.0025 


0.372 


0.009 


0.24 


0.398 


0.003 


0.384 


0.010 


0.3 


0.408 


0.003 


0.399 


0.009 


0.35 


0.421 


0.0035 


0.410 


0.010 


0.40 


0.433 


0.004 


0.428 


0.010 
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FIG. 1. An examples of J{4>) versus 4> m the symmetric sector (a) and for the broken symmetry 
sector (b). The stars (*) correspond to the values of <p at J = 0. For these results TV = 20. 
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FIG. 2. An example of J(4>) versus <p in the broken sector using the constraint effective potential. 
For these results N = 20. 
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Fig. 3 
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FIG. 3. The plot of A r versus A in the symmetric sector using lattice perturbation theory (solid 
line) and Eq. (40) (stars) and the VSM (diamonds) with m 2 = 0.1 and N = 20. 
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FIG. 4. The plot of A r versus A in the broken symmetry sector A using lattice perturbation theory 
(solid line) and VSM (diamonds) with fh 2 = —0.1 and N = 20. 
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Figure 5 
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FIG. 5. The plot of G = X r /rh r versus [3 = A + 100/A with strong coupling expansion results 
(solid line), using Eq. (40) (stars) and the VSM results (diamonds) with m% = 0.078 ±4% and N = 20. 



1 5 



10 







Figure 6 
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FIG. 6. The plot of G = \ r /m r versus (3 = (A + 100)/A with strong coupling expansion results 
(solid line), using Eq. (40) (stars) and the CPU method results (diamonds) with m%. = 0.078 ± 4% 
and N = 20. 
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